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An upper bound on the mixing efficiency is derived for a passive scalar under the influence 
of advection and diffusion with a body source. For a given stirring velocity field, the 
mixing efficiency is measured in terms of an equivalent diffusivity, which is the molecular 
diffusivity that would be required to achieve the same level of fluctuations in the scalar 
concentration in the absence of stirring, for the same source distribution. The bound on 
the equivalent diffusivity depends only on the functional "shape" of both the source and 
the advecting field. Direct numerical simulations performed for a simple advecting flow 
to test the bounds are reported. 



1. Introduction 

In this work we apply some recent developments in the analysis of the Navier-Stokes 
equations (Doering & Foias 2002) to mixing and the advection-diffusion equation. Mix- 
ing phenomena are ubiquitous with applications in atmospheric science, oceanography, 
chemical engineering, and microfluidics, to name a few. Here we focus on the generic 
problem of the advection-diffusion equation with a source that replenishes the variance 
of the passive scalar. The stirring is effected by a specified velocity field, which may or 
may not be turbulent. Our analysis of an idealised model lends mathematical precision 
and rigour to conventional scaling arguments often invoked for these kinds of problems. 

For the passive scalar, complicated behaviour — and efficient mixing — is often observed 
even for laminar velocity fields. This is the well-known effect of chaotic advection (Aref 
1984; Ottino 1989). Thus we can choose the stirring (the advecting velocity field) to 
be any divergence-free, possibly time-dependent flow field. The mixing efficiency then 
depends on specific properties of the stirring field as well as the manner in which the 
scalar concentration is injected, which is exactly what would be expected. The bound 
on mixing efficiency derived in this paper has that feature: it depends on the stirring 
field and the source distribution. This is very helpful as it allows for comparison of the 
relative effectiveness of various stirring scenarios for, say, a specified source. The bounds 
we obtain are also valid for turbulent fiows, as we make no assumptions the smoothness 
of the stirring field. A recent study by Schumacher, Sreenivasan & Yeung (2003) has 
produced bounds on the derivative moments of the concentration field; here we shall 
focus on the undifferentiated quadratic moment. As will become evident, the methods of 
this paper can also be extended to produce bounds on derivatives of the concentration 
field. 
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2. System Description 

The advection-diffusion equation for the concentration 9(x,t) of a passive scalar is 

dtO + u- Ve = kM + s , (2.1) 

where k, is the molecular diffusivity and s{x,t) is a source function with zero spatial mean. 
The domain we consider is a periodic box of side L, i.e., x G T"^, the d-dimensional torus. 
The velocity field u{x,t) could be obtained by solving Navicr Stokes or some other set 
of equations, but here we shall simply consider it to be an arbitrary divergence-free 
vector field. Hence without loss of generality we may take the solution 6{x,t) to be 
spatially mean zero at all times. 

Variations in the source term in (2.1) maintain the inhomogeneity of the concentration 
field. The stirring term may lead to the formation of sharp gradients of concentration 
that then enhance the effect of molecular diffusion. For dcfinitencss we assume that both 
the source and the stirring act on a comparable scale, £ ^ L. Because of periodicity, L/£ 
is an integer. We introduce these two distinct scales in order to be able to consider the 
infinite volume limit, L ^ oo at fixed £, for the final results. 

We shall use the fluctuations in the concentration as a useful measure of the degree of 
well-mixedness, as has long been the practise (e.g. Danckwerts 1952; Edwards, Sherman 
& Breidenthal 1985; Rehab, Antonia, Djenidi & Mi 2000). To characterise the fluctuations 
in 9, we use the variance, 

e^:=(^L-^\\9\\l,^^,^^), (2.2) 

of the spatially mean-zero concentration. The angle brackets (•) denote a long-time av- 
erage, which we will assume exists for the quantities of interest, and |H|£,2(Td) is the L'^ 
norm on T''. As control parameters we use the variance of the source and a measure of 
the kinetic energy density of the stirring field, 

:= {l-'^ , := . (2.3) 

Thus, O, S, and U arc spatio-temporal averages respectively of fluctuations in the scalar 
concentration 9{x,t), the source s{x,t), and the fluid velocity u{x,t). An efficient mixing 
configuration would have small © for a given S and U, indicating a steady-state with 
small variations in the concentration. In general we expect that increasing U at fixed S 
should decrease 8, for this represents more vigorous stirring, while increasing S at fixed U 
should augment 6. We will show in this paper that G has a lower bound proportional 
to S£/U, so that a source with large fluctuations necessarily produces a poorly mixed 
state unless U is increased sufficiently. 

In order to keep track of the effects of the amplitudes of the source variation and 
stirring intensity and their characteristic length scales independently from the influence 
of the particular "shapes" of the input and mixing functions, we decompose s and u 
into the dimensional amplitudes {S and U) and dimensionless shape functions ($ and 
T) according to 

six,t) = S^x/£,t/T), (l-'^||$||^.(t^)) = 1, (2.4) 

and 

u{x,t) = Ur{x/£,t/T), (l-'^||T||^.(t^)) = 1, (2.5) 

where r is an appropriate time scale characterising the source and stirring. Of course 
either or both may be time-independent, but in any case we presume periodicity or 
statistical stationarity with identifiable periods or relaxation times. 
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3. The Bounds 

Now consider an arbitrary smooth (dimensionless) spatially periodic function \l/(a;/£, t/r) 
normalised such that 



S = - 



*(a;/£, t/r) ^x/i, t/r) d^x ) = 1, (3.1) 

For example because of the normalisation in (2.4), ^E" = $ could be a possible choice if it 
is sufficiently smooth. Multiply (2.1) by \1/ and space-time average. Using (2.4) and (3.1) 
and integrating by parts, we may express S as 

L-'^ I (5t* + tx- V* + kA*) 6'd''a;\. (3.2) 

Note that the operator acting on ^' in (3.2) is the adjoint of the advection-diffusion 
operator, which suggests how the method can be generalised to other linear operators 
with a body source (e.g., the magnetic induction operator of dynamo theory (Childress 
& Gilbert 1995)). 

The Cauchy-Schwartz inequality implies the bound 

S < (l-'^ ll^t* + u • V* + K A*||^2(T.))'^^ e. (3.3) 

Then substituting the scaled variables T = tjr and y = x/i and using (2.5), we have 

Ue I , \i/2 

(ii^^irL^d..)) (3.4) 

where I = [0, 1] is the unit torus and 

0(y, T) := -St dr^y, T) - T(y, T) ■ V«*(y, T) + ^ {-Ay^iv, T)). (3.5) 

Here the Peclet number is Pe = U£/k. If the velocity field is time-dependent with 
timescale r, the dimensionless number Sr := £/Ut may be regarded as a Strouhal num- 
ber; in any case, we shall refer to it as the Strouhal number even if the timescale r is 
unrelated to u. 

In principle inequality (3.4) could be sharpened by varying ^ to provide as tight a 
bound as possible, as performed by Doering, Eckhart & Schumacher (2003) for the power 
consumption rate in the Navier-Stokes equations. We will not pursue that direction here; 
rather we will produce explicit limits via simple estimates. 

Applying the Minkowski inequality to (3.4), we see that 

S < ^ (ci + Pe-' C2) (3.6) 



„2 



where 

ci := (^WSrdr^ + T • Vy^.^jd)) ' , (3.7a) 

C2 := {\\AynU,,^y^\ (3.7&) 

are dimensionless constants, independent of Pe and 6. The constant Ci depends on 
dimensional quantities only through the Strouhal number; it also depends explicitly on 
the stirring shape- function T. Note also that the function ^! depends indirectly on the 
source shape-function $ through its normalisation (3.1), so that both the source and 
stirring shapes enter the bound. The constant C2 controls the diffusive part while k, only 
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enters through the Peclet number in (3.6). We still have the freedom to choose \i/ to 
optimise Ci for a particular problem, that is, for particular source and stirring shapes $ 
and T. 

For small Pe, we can focus on the C2 term in (3.6) and obtain the bound S ^ C2Qk/£'^. 
As we increase the source amplitude S, holding the other parameters constant, the time- 
averaged variance must eventually increase. An increase in the variance implies that 
the scalar is more poorly mixed. There is no avoiding this unless we increase k or decrease 
the scale of the source £: the efficiency of mixing is intrinsically related to the diffusive 
mixing rate on the scale of the source variance injection, i.e., nji"^. 

For large Pe, the more interesting limit for many physical problems, we focus on the c\ 
term in (3.6) to get the bound S < c\UQ/£. (This is true for sufficiently smooth ^ .) As we 
increase the source amplitude S, holding everything else constant, the bound (3.6) again 
implies we must eventually see an increase in the steady-state variance, 0^. However, 
unlike the small Pe case, we can now (potentially) postpone that increase by raising U, 
i.e., by stirring more vigorously. The exact value of ci depends on both shape-functions, 
but (3.7a), where ci is defined, can be broken up by the Minkowski and Holder inequalities 
to give 

ci ^ Sr (^\dTnl-(i^))^'^ + sup iVj;*! , (3.8) 

which is uniform in the shape of the stirring. The large Pe bound has the nice fea- 
ture of being independent of the diffusivity k, a result expected to hold for the passive 
scalar under turbulent or chaotic mixing. However, the linear scaling with U in Eq. (3.6) 
is not always appropriate, as will be seen in Section 5 for the specific case we have 
studied numerically. Note also that the bound (3.8) on ci still involves the velocity for 
time-dependent through the Strouhal number. If it is possible to choose ^' to be 
time-independent and still satisfy the normalisation condition (3.1) — for example if the 
source s is time-independent — then we have the bound 

ci ^ sup|Vj;*| , (3.9) 
y 

which is satisfied for all possible stirring flows {i.e., any shape function T) independently 
of U. 

We can also derive a lower bound for S. The average variance dissipation rate, e, 
satisfies 

e = Ik L-'' II V^||i.(T<*)) = (l-'' j^^ s{; x) 6{; x) d'^x'^ (3.10) 

where we have used the fact that ||^||i2(Tti) is uniformly bounded in time, which is true 
under the physical assumption that ||s||j;,2(T<i) is itself uniformly bounded in time. By 

using Poincarc's inequality in (3.10) we have (c/k)^/^ ^ (27r8)/L, and the Cauchy- 
Schwartz inequality along with the normalisation of $ in (2.4) gives SQ. Together 
these give the bound 

S ^ {2n/Lf kQ. (3.11) 

This lower bound reflects that no matter how wc stir or if wc do not stir there is still 
some diffusive dissipation of the scalar variance. The lower bound (3.11) also implies 
that if there is any variance 0^ present at the steady state, then it must be due to some 
minimum amount of amplitude of the source; stirring alone can never generate scalar 
variance in this kind of model. 

The consequence of the two bounds for S is that larger 6 must eventually imply large S 



Bound on Mixing Efficiency 5 

(from (3.11), at fixed k and L), but large S does not necessarily imply large 0, as the 
difference can be made up by a large U in (3.6). This is what makes enhanced mixing 
possible. 

We may also estimate the typical size of small scales in the scalar field. Using the 
bound S ^ e/9 mentioned above, we can transform (3.4) and (3.6) into upper bounds 
for e, viz. 

KB^^e^^ic+Pe-'c^), (3.12) 

where the lower bound is obtained via Poincare's inequality. If we define a scalar dissi- 
pation scale A, 

._ (IIV^Ili^(T^)) e 

(the Batchclor scale [Batchelor 1959], an analog of the Taylor microscale for the Navier- 
Stokes equations) then 

L/(27r) ^ A ^ £(ciPe + C2)~^/^ (3.14) 

For large Pe, the smallest possible size of this dissipation scale is proportional to Pe~^^'^, 
a standard theoretical estimate (Childress & Gilbert 1995). 



4. Mixing Efficiency and Equivalent Diffusivity 

As a physically meaningful measure of mixing efficiency, we define the equivalent dif- 
fusivity] 

Keq := /3 -Q- Ci C/^ + C2 K, (4.1) 

The factor /3 is the norm of the solution of the purely diffusive problem. 



/?:= 



(5rPe5T + A^)-i$ (4.2) 

L2(I<*) 



and the constants c\ and ci arc respectively c\ and ci multiplied by 3. The extra factor 
of /? ensures that Kcq = k for [/ = 0, which is the purely diffusiv(^ case. This corresponds 
to the choice * = {Sr Pe dr + Ay)"^^/^^ foj. ^hich C2 = 1. Note that {Sr Pe dr+Z^yY^ 
is defined in the Galerkin sense on the Fourier expansion of 

The equivalent diffusivity Koq compares the source amplitude {S) to the steady-state 
fiuctuations in the concentration field (0); as its name implies, it may be regarded as 
the molecular diffusivity needed to give a comparable amount of mixing in the absence 
of fiow. A high Peclet number steady-state mixing device should operate with as high an 
equivalent diffusivity as possible compared to the molecular diffusivity. Alternatively, we 
may interpret the ratio kIJk^ as the supression factor for the solution's variance. That 
is, if 9q is the solution of the diffusion equation with the same source but no stirring and 
Oq is it's variance, then the definition (4.1) is simply nljn"^ — 6g/0^. 

In the regime of small U, the variance is proportional to the amplitude of the source, a 

f Wc refrain from calling Kcq an "effective" diffusivity because this already carries a defini- 
tion in tlie literature (e.g. Young 1999). Tiiere tiie effective diffusivity is defined in terms of 
a large-scale gradient in the concentration, wliereas here we use the amplitude of the source, 
which makes more sense in the present context. The relationship between that traditional effec- 
tive diffusivity Kcft and Koq is Kcq = i^cti {Q/G€f , where G is a linear gradient of concentration 
(Schumacher, Sreenivasan & Yeung 2003). Other notions of effective diffusivity are also used in 
the context of anomalous diffusion (e.g. Isichenko 1992) and turbulence (e.g. Pope 2000). 
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(a) (b) 

Figure 1. The sine flow (5.1) for (a) the first half and (b) the second half of each period, 
represented here with random phases xi and X2- The velocity field alternates direction, but the 
source distribution (as represented by the shaded background) is stationary. 

response we expect when the stirring does not play an important role. A large equivalent 
diffusivity means that we are getting a well- mixed distribution (small Q) compared to 
the initial inhomogeneity in the source {S); as explained in Section 3, for fixed k and £ 
this can only be achieved by increasing U. 

The equivalent diffusivity can also be bounded from below by using (3.11), 

K^q^ KP{2we/Lf . (4.3) 

The worst lower bound for the mixing efficiency would be achieved by injecting scalar 
variance at scale £ while stirring to keep the dominant scale of the concentration fluctu- 
ation field as L. 



5. Bounds for the Sine Flow 

As an example application, we consider the well-studied two-dimensional Zeldovich 
sine flow, or random wave flow (Pierrehumbcrt 1994: Antonsen et al. 1996). This flow 
consists of alternating horizontal and vertical sine shear flows, with phase angles Xi 
and X2 G [0, 2-k] randomly chosen at each time period, r (see Figure 1). In the flrst half 
of the period, the velocity fleld is 

w(^)(a;,i) = \/2f/(0 , sin(27ra;i/L + xi)) ; (5.1a) 

and in the second half-period it is 

u^'^\x,t) = \/2f/(sin(27ra;2/i + X2) , 0). (5.16) 

The flow is incompressible, and U is deflned consistently with (2.3), so that T is read 
off from (5.1) by dropping U and replacing x/L by y. As a source function, wc choose 
s{x) = \/2S sin(27ra;i/L), from which $(y) = \/2 sin(27r2/i). Here the source and stirring 
scale length £ is equal to the system size L. The purely-diffusive solution with this source 



Bound on Mixing Efficiency 



7 




Pe 

Figure 2. Ratio of equivalent diffusivity Keq to molecular diffusivity k for the sine flow (5.1). 
The upper solid line is the upper bound (5.3) and the lower bound (4.3) is shown as a horizontal 
solid line. The upper limit (3.9) indicated by the dotted line is valid for any stirring flow with 
this source distribution. The dashed line is the result of direct numerical simulations with U 
and r fixed (i.e., with constant Sr = 1). The daslied-dot curve plots simulation data with k and 
T held constant while varying U (in this case Sr = Pe~^). 



distribution gives /? = l/(27r)^ in Eq. (4.2), and hence the lower bound Keq ^ for the 
equivalent diffusivity. 

The challenge now lies in choosing ^ to optimise the bound as best we can. The 

simplest choice is to take '5 = $, as this autornatically satisfies the norinalisation (3.1). 
Inserting that form into (3.7) (with dr'if = 0), we find Ci = v^tt and C2 = (27r)^, for a 
bound on the equivalent diffusivity 

Pe 

€—^ + 1. (5.2) 



K 



2\/27r 




We can get a tighter bound by using (3.4), which doesn't use the Minkowski inequality, 
and exploiting the statistical isotropy and homogeneity of the flow: 



(5.3) 



The bound (5.3) is actually optimal over time-independent \E' for our choice of stirring 
and source shape functions. Because T is discontinuous in time (which is not an obstacle 
to the bounding procedure), this particular velocity field does not yield a form of CI that 
is easily optimized over time-dependent So for this example with a steady source and 
a time independent multiplier ^, our bound is uniform in the Strouhal number Sr. 

Figure 2 shows the upper and lower bounds together with the results of direct numerical 
simulations of the advection-diffusion equation (2.1) with this single- mode source and 
the sine flow (5.1). The upper and lower bounds (5.3) and (4.3) are plotted as solid lines. 

There arc various ways of varying Pe in this model. For one, we can hold U and r 
fixed (both at the value 1) and vary k, in which case the Strouhal number is fixed. The 
dashed line in Figure 2 is the result of the simulation with Sr = 1 (in all simulations 
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reported here, £ = L = 1). We see that the bound quahtatively captures the behaviour 
of the equivalent difFusivity in this case, although we do not have a wide enough high Pe 
range to dctcrniine if the high-Pe asymptotic scaling agrees with the bounds'. 

Another simulation strategy is to hold k and r fixed (in this case at 10~^ and 1, 
respectively) and vary U. The data from this simulation with Sr = Pe~^ are plotted as 
the dash-dot line in Figure 2. We observe that this method of stirring is less effective 
at suppressing variance in the concentration at any value of Pe ^ 10^ where all the 
parameters coincide. In particular, with this stirring the enhancement of the equivalent 
diffusivity tends to saturate rather than increase indefinitely at high Pe. This can be 
understood as a "resonance" efi'ect of the periodic boundary conditions: as U increases 
at fixed T and the typical displacement Ut exceeds i, the velocity field merely maps 
the concentration fluctuations onto periodic copies of themselves rather than mixing 
it up within each periodic cell. Because the bound in (5.3) is uniform in r, it should 
be compared at a given value of Pe against the largest possible equivalent diffusivity 
achievable by any r, which can only improve the agreement with the bound. These 
different simulation schemes illustrate the importance of the Strouhal number for the 
mixing efficiency and, not unexpectedly, highlight the need for further analysis to extract 
the Strouhal number dependence of the best bounds. 

The dotted line in Figure 2 is a weaker upper bound obtained from (3.9) using \1/ = i>; 
it sets an absolute limit on the mixing efficiency achievable with any stirring field shape 
at any Strouhal number for this particular source distribiition. 

Finally, we note that there are flow flelds at arbitrarily high values of Pe and arbitrary 
Sr that saturate the lower bound Keq/K = 1 for this source shape. Indeed, any flow 
field u with no X2 dependence (and arbitrary xi and t dependence) simply moves the 
scalar along iso-concentration lines to no effect whatsoever. This simple example is a 
particular case of a more general result concerning the existence of "ineffective" stirring 
fields (Young 2004) — essentially integrable fields without chaos. 

6. Conclusions 

It is encouraging that the equivalent diffusivities in Figure 2 rise away from the diffusive 
lower bound as Pe increases, indicating that there is hope of more nearly saturating the 
upper bound with more complex flows. From the example of the sine flow, it is clear that 
in general there is a nontrivial Sr dependence, even for steady sources, that warrants 
further investigation. Unlike the solution of the full problem which requires a nonzero 
diffusivity to keep O uniformly bounded in time, the bounding procedure does not require 
any diffusivity. That is, for large Pe we may neglect C2 from the bound altogether and 
focus on ci to try and minimize it with respect to Of course, the resulting optimal 
bound on Kcq/U£ may still depend in a complicated way on Pe and Sr for specific stirring 
and source distributions. 

The high-Pe scaling of the bound obtained in this paper might be related to an anal- 
ogous one in combustion theory (Constantin et al. 2000). There it was fomid that the 
bulk burning rate V can satisfy an "optimal linear enhancement bound," V ^ KU , 
where K is & constant and U is the magnitude of the advecting field. The type of flow 
required for linear enhancement, called "percolating flows" in Constantin et al. (2000), 
connects distant regions of unburned material. Perhaps these flows also produce linear 
asymptotic scaling with Pe for the equivalent diffusivity enhancement, but we have not 
yet investigated this. 

Although we specifled a body source in our problem with periodic conditions, a source 
of concentration at impenetrable boundaries can be mimicked by a sharp source concen- 
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trated near the walls (Balmforth & Young 2003). However, the type of wall boundary 
condition that can be modelled in this manner is restricted to fixed scalar flux. 

In closing wc note that all of our analysis, as well as the general result that Kcq/ k := 
PSf' / hiQ ^ ciPe + C2, depend on the source distribution being smooth enough to have 
a finite variance 5^. Point sources, for example where s ^ 5{x)^ may be of interest in 
applications but do not have finite variance. In this situation we may still define the 
mixing efficiency and an equivalent diffusivity via Kcq/n := Oo/0 where 0^ and Oq are 
the scalar variances with and without the stirring; these scalar variances are finite even for 
^-like sources in two and three spatial dimensions. However the anticipated behaviour 
suggested by the consideration of smooth sources, i.e., that the equivalent diffusivity 
enhancement Keq/K and/or its upper bound could be ^ Pe, may not be realized with 
more singular sources. The investigation of those models is left for future work. 

We thank P. Constantin, K. R. Sreenivasan, and W. R. Young for helpful comments. 
J.-L.T. and C.R.D. are grateful for the hospitality of the 2002 Summer Program in 
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